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Abstract 

The paper is concerned with sharp estimates of constants in classical Poincare inequalities and Poincare-type 
inequalities for functions having zero mean value in a simplicial domain or on a part of its boundary. These esti¬ 
mates are important for quantitative analysis of problems generated by differential equations, where numerical 
approximations are typically constructed with the help of simplicial meshes. We suggest easily computable rela¬ 
tions that provide sharp bounds of the respective constants and compare these results with analytical estimates 
(if they are known). In the last section, we present an example that shows possible applications of the results 
and derive a computable majorant of the difference between the exact solution of a boundary value problem 
and an arbitrary finite dimensional approximation computed on a simplicial mesh, which uses above mentioned 
constants. 


1 Introduction 

Let T be a bounded domain in {d > 2) with Lipschitz boundary dT. It is well known that the Poincare inequality 
([321133]) 

IIv^IIt ^ ||Vw||t' (1) 

holds for any 

W e H^T) := e H\T) \ {1^1}.^ = 0 }> 

where ||w||t denotes the norm in L^(T), {|w|}y := fj,wdx is the mean value of w over T, and Rl is the Lebesgue 

measure of T. The constant depends only on T and d. 

Poincare-type inequalities also hold for 

w e i?yr,r) := {a; G H^iT) \ {|n;|}r = 0 }, 

where P is a measurable part of dT such that meaSd-iP > 0 (in particular, P may coincide with the whole boundary). 
For any w G iT^(T, P), we have two inequalities similar to (|T|). The first one 

Ikllr <DP||Vn;||T (2) 

is another form of the Poincare inequality o, which is stated for a different set of functions and contains a different 
constant, i.e. < Cp. The constant Cp is associated with the minimal positive eigenvalue of the problem 

—Am = Am in T; 9„m = A{|m|}p on P; 9 „m = 0 on i9r\r; VuGff^(T,r). (3) 

We note that inequalities of this type arose in finite element analysis many years ago (see, e.g., 0), where @ was 
considered for simplexes in R^. The second inequality 

Ikllr < UfJlVwlR (4) 
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estimates the trace of w S H^{T,r) on F. It is associated with the minimal nonzero eigenvalue of the problem 

-Am = 0 in T; dnU = Xu on T; dnU = 0 on dT\T-, yu€H\T,r). (5) 


The problem ([5]) is a special case of the Steklov problem [39] , where the spectral parameter appears in the boundary 
condition. Sometimes ([S]) is associated with the so-called sloshing problem, which describes oscillations of a fluid in 
a container. Eigenvalues and eigenfunctions of the sloshing problem have been studied in [UlllITllllllKll] and 
some other papers cited therein. 

Exact values of Cp , Cp"^, and Cp are important from both analytical and computational points of view. Poincare- 
type inequalities are often used in analysis of nonconforming approximations (e.g., discontinuous Galerkin or mortar 
methods), domain decomposition methods (see, e.g., [El [12] and [41]), a posteriori estimates [37], and other ap¬ 
plications related to quantitative analysis of partial differential equations. Analysis of interpolation constants and 
their estimates for piecewise constant and linear interpolations over triangular finite elements can be found in [23] 
and literature cited therein. Einally, we note that [5] introduces a method of computing lower bounds for the eigen¬ 
values of the Laplace operator based on nonconforming (Crouzeix-Raviart) approximations. This method yields 
guaranteed upper bounds of the constant in the Eriedrichs’ inequality. 

It is known (see [30] 1 that for convex domains 


^ diam(r) 

^ — TT 

Eor triangles this estimate was improved in [22] to 


^ diam(T) 
— ii,i ’ 


where ji_i ~ 3.8317 is the smallest positive root of the Bessel function Ji. Moreover, for isosceles triangles from 
it follows that 


Ct < := diam(T) • 


tan(f))” 


a e (0, |], 
a e (|, f], 

a e (|,7r). 


( 6 ) 


Here, jo,i ~ 2.4048 is the smallest positive root of the Bessel function Jq. A lower bound of for convex domains 
in was derived in [E] and it reads 

CP > diam(T) _ .j. 

^ — 2jo,i ^ ' 


Analogously, work [2T] provides lower bound 


> 


P 

4 TT ’ 


( 8 ) 


which improves 0 for some cases. Here, P is perimeter of T. 

In [22], exact values of Cf and are found for parallelepipeds, rectangles, and right triangles. Subsequently, 
we exploit the following two results: 


1. If T is based on vertexes A = (0,0), B = {h,0), C = (0,h) and T := {xi G [0,/i], X 2 = O} (i.e., T coincides 
with one of the legs of the isosceles right triangle), then 


Cf = A and C;?’' 


( . . V 

VCo tanh(Co)/ 


where Co and Co are unique roots of the equations 


zcoi{z) -1-1 = 0 and tan(z) -I- tanh(z) = 0, 


(9) 


( 10 ) 


respectively, in the interval (0,7r) . 

2. If T is based on vertexes A = (0,0), B = (/i,0), C = (|, |), and T coincides with the hypotenuse of the 
isosceles right triangle, then 

CF = 4 and C?■=(|)'/^ 


It is worth emphasizing that values of for right isosceles triangles follow from the exact solutions of the Steklov 
problem related to the square. This specific case was discussed in the work [M]- 

Exact value of constants in the classical Poincare inequality are also known for certain triangles: 
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1. For the equilateral triangle based on vertexes A = (0,0), B = (1,0), C = ( 5 ,^), where 

r := {xi G [ 0 , 1 ]; X 2 = 0 }, the constant 

r'P — JL 

f, Va 47r 

is derived in m- 

2 . For the right isosceles triangles based on vertexes A = (0,0), B = (1,0), C = (i, i) and based on 
A = (0,0), B = (1,0), C = (0,1), we have 


_ 1 

^f,V4 V2ir 


and 



TT ’ 


respectively. Proofs can be found in m and |28) . 

Explicit formulas of the same constants for certain three-dimensional domains are presented in papers m and [TS] . 

The above mentioned results form a basis for deriving sharp bounds of the constants Cp, C'p'', and for 
arbitrary non-degenerate triangles and tetrahedrons, which are typical objects in various discretization methods. In 
Section [21 we deduce guaranteed and easily computable bounds of Cp, CjF^, and for triangular domains. The 
efficiency of these bounds is tested in Section [31 where Cp and Cj? are compared with lower bounds computed 
numerically by solving a generalized eigenvalue problem generated by Rayleigh quotients discretized over sufficiently 
representative sets of trial functions. In the same section, we make a similar comparison of numerical lower bounds 
related to the constant with obtained upper bounds and existing estimates known from [2T1 [22] and [TOj . Lower 
bounds of the constants presented in Section [3] have been computed by two independent codes: the first code is 
based on the MATLAB Symbolic Math Toolbox [JD], and the second one uses The FEniCS Project [21]. Section 
ID is devoted to tetrahedrons. We combine numerical and theoretical estimates in order to derive two-sided bounds 
of the constants. Finally, in Section [5] we present an example that shows one possible application of the estimates 
considered in previous sections. Here, the constants are used in order to deduce a guaranteed and fully computable 
upper bound of the distance between the exact solution of an elliptic boundary value problem and an arbitrary 
function (approximation) in the respective energy space. 


2 Majorants of C'p and Cp’^ for triangular domains 

Let T be based on vertexes A = (0, 0), H = {h, 0), and C = [hp cos a, hp sin a) and 

r:={a:iG[0,/i]; X2 = 0}, (11) 

where p > 0, /i > 0, and a £ ( 0 , 7 r) are geometrical parameters that fully define a triangle T (see Fig. [T]). Easily 
computable bounds of Cf and C'p'' are presented in Lemma [T] below, which uses mappings of reference triangles to 
T and well-known integral transformations (see, e.g., [TTjL 



Figure 1: Simplex in R^. 


Lemma 1 For any w G iL^(T, F), the upper bounds of constants in the inequalities 


|w||t < Cp/i||Vr(;||T and ||w||r < Vu>||t 


( 12 ) 


are defined as 

->? —P 


Cf < = min {^f,, , 7.*}. | «n<i Cp" < Cr" = min ( 7 ^', Cp",., 7?/. G 

respectively. Here, 


;Tr 


lin 


rr I 
gAJ’ 


7v, = , 7v, = Mvl, 7v. = {P sin «) '" 7?/,, 7^4 = (2p sin a) ' ^ 7 ?/^, 


“V 2 


,-V 2 p 
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where 


{p, a) =\(^ + + {l + + 2 cos 2a) j, 

1 j 

p^/^{p, a) = 2p^ — 2p cos a + 1 + ((2p^ + l)(2p^ + 1 — 4p cos a + cos 2a)) 

and CpP , ss 0.49291, ^ 0.65602 and , « 0.24646, « 0.70711, where f is defined as 

r.-/2 ’ r,V2 r-/4 ’ r.V4 ’ ■' 

f := {xi e [0,1]; X2 = o}. 

Proof: Consider the linear mapping ^ T with 

h ph cos a 


X = (a;) = i, where Bi^j^ = 

For any w e we have the estimate 


0 ph sin a I ’ 


deti?7r/2 = ph'^ sin a. 


w 


<Cf^, II Vw 


■T,/, - -r-/ 




(13) 

(14) 

(15) 


(16) 


where is the constant associated with the basic simplex T-Kj^ based on 4 = (0,0), B = (1,0), and C = (0,1), 

Note that 


'^1'2 


1 i|zc |'2 


ph"^ sin a II ^ 117"’ 


and 


^ J A.f^{h,p,a)VwVwdx, 


(17) 

(18) 


where 


A.u(h,p, a) = h^ 


1 + cos^ a 


p^ sin a cos a 


p^ sin a cos a p^ sin^ a 


It is not difficult to see that Xma^{A7,i^) = h^p^/^{p,a), where p^i^{p,a) is defined in (fT^ . From (IT6l) . (flTl) . and 
m, it follows that 

II^IIt < II Vw||t, ^^/^{p,a) = pj^l{p,a). (19) 

Notice that w G H^{T,r) yields 


llrcllp := J w{x) ds = h J r(;(a;(i))ds = h J w)ds = 0. 


Therefore, above mapping keeps w € B^(T,r). 

In view of inequality (|3]), for any w G H^(T-,^/^,r) we have 

||*||f<Cp^VJ|Vii||f^^^, 

where is the constant associated with the reference simplex Since 

II ^ lip = ill Hr, 


we obtain 


Ir < V„.|lr, 7.7 Ja«) = 


( 20 ) 


Now, we consider the mapping ^ T, where Ti,/^ is based on A = (0,0), B = (1, 0), and C = (^, 2 ), i-®-, 


X = Tirjjyx) = i?7r/^ i, where = 


h 2phcosa-h' 
0 2phsina j 


det B„l^ = 2ph?' sin a, 
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which yields another pair of estimates for the functions in H^{T,r): 


and 


w\\t < 7-/4 C'f,V4 ^ II 


( 21 ) 


Ir < 7 . 7 , 7",>A<Ao) = . 


( 22 ) 


where is defined in (IT4l) . Now, (IT^ follows from (IT^ . (l20l) . (1^ . and (l22l) . 


□ 


Analogously to Lemma [TJ one can obtain an upper bound of the constant in ([T]). For that we consider three 
reference triangles (defined earlier), and based on vertexes A = (0, 0), B = (1,0), C = (i, -^). 

Lemma 2 For any w € F{^(T), the constant in 


IkllT < cEhw^wWt, 


(23) 


is estimated as 


Cl<cl = min {yP ^ ^ y?/, ,/, }■ 


(24) 


Here, y^/ = yP /3 = pif, where and being defined in (fT^ and (fT^ . respectively, 


and 


a) = |(1 + — p cos a) + 2 (i(l + p^ — p cosa)^ — \p^ sin^ a) 


V 2 


(25) 


^fy/4 = 


Proof: The mapping : Ttt/^ —>• T coincides with © from Lemma [H It is easy to see that w € H^{T) provides 
that w € H^{T). The estimate 


I W IIt < yw, CS h II Vw IIt, X^/Sp,a) = pJ/liP^a) 


V2 '-"f 7/2 II ll'P’ '^V 2 yf'’ - ^V 2 

is obtained by following steps of the previous proof. From analysis of mappings 

h 


(26) 


X = = Bir/g X, where = 


and 

X = •^V 4 (^) = where 

we obtain alternative estimates 


h -^( 2 pcosa — 1 )-/i' 
,0 :^^sina ^ 


, detSTT/g = ^ sin a > 0 , 


h 2ph cos a-h\ , , ^7 <772 


0 2 p/isina 


deiB^/^ = 2ph sino > 0 , 


w IIt < yP /3 h II Vw IIt, yP/ 3 (p,a) = p!fJy.P.o^) 

w^IIt < X^/C^„/h\\Vw\\T, X^/AP,a) = pjfi{p,a) 


V 4 ,v. 


V 4 


(27) 

(28) 


where p^/,^{p,a) and p7r/^(p,a) are defined in (PSI) and (|TT)) . respectively. Therefore, follows from combination 
of (PSI) . (I?7l) . and (1^ . □ 


3 Minorants of and for triangular domains 

3.1 Two-sided bounds of Cf and Cj? 

Majorants of Cf and Cj? provided by Lemma [T] should be compared with the corresponding minorants, which can 
be found by means of the Rayleigh quotients 


77.p [w] 


-JVHIt 

k-ll’^llrll'r 


and 


77 .p''[?t] = 


’"-{Nllrllr ■ 


(29) 
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Lower bounds are obtained if the quotients are minimized on finite dimensional subspaces C H^{T) formed 
by sufficiently representative collections of suitable test functions. For this purpose, we use either power or Fourier 
series and introduce the spaces 


:= span{ } and := span{ cos(7rzx) cosijrjy) }, 

where i,j = 0,... ,fV, (i,j) ^ (0,0) and 

dimFj^ = dim V 2 ^ = M(iV) := {N + if - 1. 

The corresponding constants are denoted by C_^’^ and where M indicates on number of basis functions in 

auxiliary subspace used. Since and are limit dense in H^{T), the respective minorants tend to the exact 
constants as M{N) tends to infinity. 

We note that 

inf 72.p[w] = inf -n— n = inf = fp 

wem{T) wem(T) ll’"-{l’"l}rllT weHi{T,r) 

Therefore, minimization of the first quotient in (0^ on or 
7^p“'[w], we apply similar arguments. 

Numerical results presented below are obtained with the help of two different codes based on the MATLAB 
Symbolic Math Toolbox [JD] and The FEniCS Project [53] . Table |T] demonstrates the ratios between the exact 
constants and respective approximate values (for the selected p and a). They are quite close to 1 even for relatively 
small N. Henceforth, we select = 6 or 7 in the tests discussed below. 


(30) 


yields a lower bound of Cf. For the quotient 



II 

p=l 


II 

0. 

N 

M{N) 

p,M,P ,„p 

-r /^r,V2 

—e ' r,V2 

p,M,P ,„p 

-r /^r,V4 

„M,Tr Tr 

-r 

1 

3 

0.8801 

0.9561 

0.8647 

1.0000 

2 

8 

0.9945 

0.9898 

0.9925 

1.0000 

3 

15 

0.9999 

0.9998 

0.9962 

1.0000 

4 

24 

1.0000 

0.9999 

1.0000 

1.0000 

5 

35 

1.0000 

1.0000 

1.0000 

1.0000 

6 

48 

1.0000 

1.0000 

1.0000 

1.0000 


Table 1: Ratios between approximate and reference constants with respect to increasing N. 


In Figs. |5a| and |5cl we depict for M{N) = 48 (thin red line) for T with p = 1, and a € ( 0 , 7 r). 


_p _p 

Guaranteed upper bounds C-„/ = 7 ?, Cf and C-^/ = 7 , 

_p r t:) A 

line illustrates Cp = min 


Analogously in Figs. [3al and [3bl a 


are depicted by dashed black lines. Bold blue 
red marker denotes the lower bound 


^M,Tr M[N) = 48) of the constant C'j?''. It is presented together with the upper bound Cpl^ (blue marker), 
which is defined as minimum of C-nj^ = ^'"/a ~ ^"^a represents this information in 

the digital form. 

Fig. [5a| corresponds to the case p = Notice that for a = ^ the constant Cf is known and the computed 


Tr 


lower bound (red marker) practically coincides with it (see, e.g., Fig. [2b| . Since in this case, the mapping 

is identical, the upper bound also coincides with the exact value. An analogous coincidence can be observed 
for and in Fig. |3a| In Fig. |2cl the red curve, corresponding to C^’^, coincides with the blue line of Cf 

at the point a = § (due to the fact that for this angle J' is the identical mapping and T coincides with (see 
Fig. [2dl)'). Fig, [^exposes similar results for Cp^’’^ and (C'^a)- 

Figs, m and ID demonstrate the same bounds for p = -^ and |. We see that estimates of Cp and C^ are very 
efficient. Namely, := G [1.0463,1.1300] for p = -^ and G [1.0249,1.1634] for p = |. Analogously, 


eff ■ c 


48,Tr 

b=Lr 


G [1.0363,1.3388] for p = ^ and G [1.2917,1.7643] for p = |. 


3.2 Two-sided bounds of 


The spaces and Vf can also be used for analysis of the quotient TZt [u>] = , which yields guaranteed 

lower bounds of the constant in O- The respective values are denoted by . These bounds are compared with 

1®_diam(T) 

O T’ — - 


Jl 


and Cp := max | 1 (see ©“(HI), respectively) as well as the one derived in Lemma |2| 
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Figure 2: Two-sided bounds of Cp for T with different p. 




Figure 3: Two-sided bounds of Cp'' for T with different p. 
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OL 

^48,P 
}=±r 

Cr 

^48,Tr 


^48,P 

Cr 

^48,Tr 

Cp- 

7r/18 

0.2429 

0.2657 

1.2786 

1.5386 

0.3245 

0.3486 

1.2572 

1.6971 

■K I'd 

0.2414 

0.2627 

0.9289 

1.0838 

0.3248 

0.3493 

0.9058 

1.2116 

7r/6 

0.2389 

0.2577 

0.7919 

0.8792 

0.3268 

0.3527 

0.7632 

1.0118 

27r/9 

0.2379 

0.2507 

0.7259 

0.7543 

0.3339 

0.3636 

0.6906 

0.9201 

57r/18 

0.2632 

0.2722 

0.6945 

0.7503 

0.3514 

0.3884 

0.6529 

0.9003 

Tl jZ 

0.3008 

0.3220 

0.6829 

0.8348 

0.3809 

0.4269 

0.6362 

0.8634 

77r/18 

0.3382 

0.3694 

0.6840 

0.8432 

0.4173 

0.4721 

0.6332 

0.7840 

47r/9 

0.3740 

0.4140 

0.6947 

0.7973 

0.4556 

0.5187 

0.6404 

0.7162 

7r/2 

0.4075 

0.4554 

0.7136 

0.7801 

0.4929 

0.4929 

0.6560 

0.6560 

h-K jd 

0.4382 

0.4933 

0.7409 

0.7973 

0.5280 

0.5340 

0.6797 

0.7162 

ll7r/18 

0.4660 

0.5165 

0.7779 

0.8432 

0.5600 

0.5710 

0.7125 

0.7840 

27r/3 

0.4905 

0.5361 

0.8274 

0.9118 

0.5884 

0.6037 

0.7569 

0.8634 

137r/18 

0.5115 

0.5552 

0.8948 

1.0040 

0.6129 

0.6318 

0.8175 

0.9607 

l-Kjd 

0.5289 

0.5720 

0.9898 

1.1292 

0.6332 

0.6550 

0.9033 

1.0874 

b-Kja 

0.5426 

0.5856 

1.1334 

1.3107 

0.6492 

0.6733 

1.0332 

1.2673 

87r/9 

0.5524 

0.5956 

1.3796 

1.6118 

0.6607 

0.6865 

1.2565 

1.5623 

177r/18 

0.5583 

0.6017 

1.9436 

2.2851 

0.6676 

0.6944 

1.7692 

2.2179 

Table 2: Two-sided bounds of Cf 

and Cj? for T for 

o. G (0, 

tt) and different 




(a) p = ^ (b) P = I 

Figure 4: Two-sided bounds of Cp for T with different p. 




(a)p=^ (b)p=| 


Figure 5: Two-sided bounds of for T with different p. 





















































Figure 6: Cy®, C^, and C]^ for T with a S (0,7r) and different p. 

In Figs. IBal I6b[ and I6d[ we present c!^’^ (in this case M{N) = 48) together with (blue think line), Cy® 
and for a G (0,7r), and p = and We see that (red thin line) indeed lies within the admissible 

_p _p 0 

two-sided bounds. From these figures, it is obvious that new upper bounds Crp are sharper than for T with 

p ^ 1. True values of the constant he between the bold blue and thin red lines, but closer to the red one, which 
practically shows the constant (this follows from the fact that increasing M{N) does not provide a noticeable change 
for the line, e.g., for M{N) = 63 maximal difference with respect to figure does not exceed le-8). Also, we note 
that, the lower bound C]^ (black dashed line) is quite efficient, and, moreover, asymptotically exact for a ^ tt. 

_p A 

Due to [H] and |3], we know the improved upper bound Crp’ (cf. ([6|)) for isosceles triangles. In Fig. |6cl we 

compare C_^’^ {M{N) = 48) with both upper bounds Cy (from the Lemma [2|) and Cy’ (black doted line). It 

is easy to see that Crp’ (black dashed line) is rather accurate and for a ^ 0 and a —>■ tt provide almost exact 

estimates. Cp (blues thick line) improves Cp only for some a. The lower bound (red thin line) indeed 
—P,A 

converges to Cp as T degenerates when a tends to 0. 

3.3 Shape of the minimizer 

Exact constants in © and 0 are generated by the minimal positive eigenvalues of ([3]) and 0. This section 
presents results related to the respective eigenfunctions. In order to depict all of them in a unified form, we use the 

3 

barycentric coordinates Aj G (0,1), i = 1, 2,3, X) = 1- 

i=l 
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Figs. [3 and [5] show the eigenfunctions computed for isosceles triangles with different angles a between two legs 
(zero mean condition is imposed on one of the legs). The eigenfunctions have been computed in the process of 
finding C_^’^ and . The eigenfunctions are normalized so that the maximal value is equal to 1. For a = ^, 

the exact eigenfunction associated with the smallest positive eigenvalue Ap = is known (see [29] 1: 

= cos(^) + cos . 

Here, Co is the root of the first equation in (fTOl) (see FielTdl). We can compare with the approximate eigenfunction 
Up ’ computed by minimization of 7^p[w] (this function is depicted in Fig. [7c]). 

Eigenfunctions related to the constant are presented in Fig. |8l Again, for a = ^ we know the exact 

eigenfunction 

M]? = cos(Coa;i) cosh (Co(a ^2 - h)) + cosh(Co xi) cos (Co(a ;2 - h)), 
where Co is the root of second equation in m (see Fig. I8d|) . This function minimizes the quotient 7^p''[?n] and 
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1 


1 


/ \ 48,Tr TT 

(a) lip ’ , Oi= Q 


(b) u 


48.Tr 

r ’ “ “ 3 




yields the smallest positive eigenvalue . it is easy to see that numerical approximation 0^’"^ (for 

M{N) = 48) practically coincides with the exact function. 

Typically, the eigenfunctions associated with minimal positive eigenvalues expose a continuous evolution with 
respect to a. However, this is not true for the quotient 7 ?.t[w], where the minimizer radically changes the profile. Fig. 
[5^ indicates a possibility of such a rapid change at a = ^, where the curve (related to C^) obviously becomes non¬ 
smooth. This happens because an equilateral triangle has double eigenvalue, therefore the minimizer of TZtIw] over 
V-^ changes its profile. Figs. l9aM9il show three eigenfunctions Ut® 2 > corresponding to three minimal 

eigenvalues ^t* 2 > ^"^cl Af^g. All functions are computed for isosceles triangles and are sorted in accordance with 
increasing values of the respective eigenvalues. It is easy to see that at a = ^ the first and the second eigenfunctions 
swap places. Table |3] presents the corresponding results in the digital form. 
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Figure 9: Eigenfunctions corresponding to C_^’^ with M = 48 on isosceles triangles T G with a = ^, 
^ + e in barycentric coordinates. 



^T,i 

/^48 \48 

^T,i ^T,i 

/^48 \48 

/^48 \48 

^r,i 

P = 1 

Ut 1 

“'T,3 

0.2419 17.0951 

0.2229 20.1216 

0.1353 54.6024 

0.2387 17.5463 

0.2387 17.5463 

0.1378 52.6396 

0.2537 15.5404 

0.2355 18.0309 

0.1422 49.4818 

II 

0 . 

Ut 1 

“'T,3 

0.23137 18.6804 

0.17082 34.2707 

0.1229 66.2058 

0.23671 17.8471 

0.17435 32.8970 

0.12789 61.1402 

0.24336 16.8850 

0.17642 32.1295 

0.13298 56.5493 

p=f 

Ut 1 

“'T,3 

0.34714 8.2983 

0.24485 16.6801 

0.18258 29.9981 

0.35523 7.9247 

0.24885 16.1482 

0.19084 27.4575 

0.3648 7.5143 

0.25125 15.8412 

0.19845 25.3921 


Table 3: C_^’^ and corresponding to the first three eigenfunctions in Fig. [51 


It is worth noting that for equilateral triangles two minimal eigenfunctions are known (see |26jl: 

Ml = cos - 1)^ - cos cos (^^{2xi - 1)^, 

U 2 = sin ^2^(2xi - 1)^ + cos (^ 2 : 2 ) sin ( 2 x 1 - 1)^ 
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w|=l 





























M{N) 


^Tr,M 

f 

r,V2,« 

QTr,M 

r , 2, a 

r,V2iQ' 

^Tr,M 

r,V2,« 

fjP,M 

^Tr,M 

7 

0.32431 

0.760099 

0.325985 

0.654654 

0.360532 

0.654654 

0.4152099 

0.686161 

26 

0.338539 

0.829445 

0.340267 

0.761278 

0.373669 

0.751615 

0.4274757 

0.863324 

63 

0.341122 

0.831325 

0.342556 

0.762901 

0.375590 

0.751994 

0.4286444 

0.864595 

124 

0.341147 

0.831335 

0.342589 

0.762905 

0.375603 

0.751999 

0.4286652 

0.864630 

215 

0.341147 

0.831335 

0.342589 

0.762905 

0.375603 

0.751999 

0.4286652 

0.864630 


Table 4: . and with respect to M{N) for Tg ^ with p = 1, 9 = ^, and different a. 

i- • / 2 5 ^ -^5 / 2 1 ^ ’ 


These functions practically coincide with the functions u^i and m |?2 presented in Fig. I9dl Finally, we note that 

this phenomenon (change of the minimal eigenfunction) does not appear for p = ^ or p = |. The eigenvalues as 
well as the constants corresponding to the eigenfunctions presented in Fig. |9] are shown in the Tabled 

4 Two-sided bounds of and C'p'^ for tetrahedrons 


X2 



X2 



Figure 10: Simplex in R^. 


Figure 11: Coordinate of the vertex D. 


We orient the coordinates as it is shown in Fig. [TU] and define a non-degenerate simplex in R^ with vertexes 
A = (0,0,0), B = (/ii,0,0), C = {0,0,hs), and D = (^2 sind cosa,h 2 sin0 sina,/i 2 cosa), where hi and /13 are 
the scaling parameters along axis and Ox^, respectively, AD = /i 2 , a is a polar angle, and 9 is an azimuthal 
angle (see Fig.[TT]). Let F be defined by vertexes A, B, and C. 

To the best of our knowledge, exact values of constants in Poincare-type inequalities for simplexes in R^ are 
unknown. Therefore, we first consider four basic (reference) tetrahedrons with /i 2 = 1, 0 = ^, and di = ^, 0:2 = f) 
da = and ai = The respective constants are found numerically with high accuracy (see Tabled which shows 
convergence of the constants with respect to increasing M{N)). Henceforth, Tg ^ denotes a reference tetrahedron, 

where 9 and d are certain fixed angles. By Tg ^ we denote the respective mapping Tg ^ : Tg ^ T. Then, for an 
arbitrary tetrahedron T, we have 

IkllT < CF/i 2 ||Vu||t and |k||r < C?-|| Vu||t (31) 


13 


























with approximate bounds 


Cf ^ Ct? = min 

<i={V4.V3.-/2dV3} 




} 


and 


C? < = 


i={V4.V3.V2dV3} 


|_^V 2 ,a ^r,V 2 ,a 


(32) 


(33) 


where . and C'^'', . are the constants related to four reference tetrahedrons from Table IH and - and 


r,”/ 2 ,a 


r.V 2 ,a 


(i (IMll l are generated by the mapping T^/ 2 ,& Here, the reference tetrahedrons are defined 

based on A = (0,0,0), B = (1,0,0), C = (0,0,1), D = (cosd,sinQ;, 0) with a = {j, ^}, and is 

presented by the relation 

X ^7r^2,a(^) -H7r^2,a^j — 1,2,3 ^2 


/ hi. 
/12 

0 
\ 0 


sin a 

sin Q sin 6 
sin OL 
cos 9 
sin a. 


where v{p, a) = cos a sin 0 — ^ cosd, det = hi h 2 in a" ^ analogy with the two-dimensional case (see 

dun), 7 )^^ ^ and 7 )iy^ ^ depend on the maximum eigenvalue of the matrix 



/ bh + 6?2 

^12^22 

& 12&32 


^12^22 

bl 2 

622^32 


V ^ 12^32 

^ 22^32 

bla + bh 


The maximal eigenvalue of the matrix is defined by the relation Amax( 447 r/ 2 _Q) = hf Pa,e,a with 

Pc.,e,&= (£5" -£38^^" + ^£i), 

where 


£1 — bh + 6^2 + ^22 + ^32 + ^331 

£2 = bh 5^2 + bfi bl^ + bj^ bj^ + bl^ bj^ + bl^ bj^, 


^3 = f - (^)^ 

£4 = {^f-^ + ^blibl,bl 
£5 = £4 + {£! + £!f^^- 


Therefore, ^ 7 ^/ 2 ,a ^ad T7ry'2^a are as follows: 


- 

1 ^/ 2 ,a 


„V 2 

f^Tr/2,ai 


- 

' l2,a 


] 

( sin 6c \ 
p sin oc sin 9 J 


P 

T’f /2 ,a ■ 


(34) 


Lower bounds of the constants Cp and Cj? are computed by minimization of 7^p[r(;] and over the set 

Vg^ C H^(T), where 

■= [y^^3k = x'^y^z'', i,j,k = 0,...,N, (z, j, fc) 7 ^ (0, 0,0) | 
and dimVi(^ = M{N) := {N + 1)3 _ 

The respective results are presented in Tables [5] and [6] for T with hi = 1, hs = 1, and p = 1. We note that exact 
values of constants are probably closer to the numbers presented in left-hand side columns. For 9 = Y 2 , we also 
present estimates of C_^'^ and (red lines) graphically in Fig. [121 


5 Example 

Constants in the Friedrichs’, Poincare, and other functional inequalities arise in various problems of numerical 
analysis, where we need to know values of the respective constants associated with particular domains. Constants 
in projection type estimates arise in a priori analysis (see, e.g., [giniEi]). Constants in Clement’s interpolation 
inequalities are important for residual type a posteriori estimates (see, e.g., [l]|42], and [7], where these constants 
have been evaluated). Concerning constants in the trace inequalities associated with polygonal domain, we mention 
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OL = 


CL = 


a = 


3 


a = 


TT 

6 


TT 

4 


TT 

2 


6 


qM,P Qp 


qM,F gp 


^M,P 


CpP 


^M,P gp 


7r/6 

0.23883 

0.49035 

0.24621 

0.49841 

0.25870 

0.51054 

0.29484 

0.51308 

7r/4 

0.23883 

0.45388 

0.24621 

0.46173 

0.25870 

0.47683 

0.29484 

0.49075 

7r/3 

0.29666 

0.41958 

0.31194 

0.42259 

0.33489 

0.43724 

0.38976 

0.46002 

7r/2 

0.34302 

0.35667 

0.34112 

0.34115 

0.34256 

0.34259 

0.37559 

0.37560 

27r/3 

0.40428 

0.41958 

0.40562 

0.42259 

0.40927 

0.43724 

0.42867 

0.46002 

37r/4 

0.42890 

0.45388 

0.43110 

0.46173 

0.43505 

0.47683 

0.45017 

0.49075 

5it/6 

0.44964 

0.49035 

0.45193 

0.49841 

0.45539 

0.51054 

0.46607 

0.51308 





27r 






OL — 

2 

OL — 

3 

OL — 

4 

OL — 

6 

0 

Cv'^ 

CP 

^M,P 

Llp 

CP 

C(^’P 

CP 

^M,P 

Llp 

CP 

7r/6 

0.29484 

0.51308 

0.33069 

0.51792 

0.34468 

0.52253 

0.35499 

0.52694 

7r/4 

0.29484 

0.49075 

0.33069 

0.50261 

0.34468 

0.51308 

0.35499 

0.52253 

TvfS 

0.38976 

0.46002 

0.43880 

0.48413 

0.45742 

0.50261 

0.47106 

0.51792 

7^/2 

0.37559 

0.37560 

0.42865 

0.42867 

0.45017 

0.45731 

0.46607 

0.47811 

27r/3 

0.42867 

0.46002 

0.45997 

0.48413 

0.47457 

0.50261 

0.48598 

0.51792 

Znli 

0.45017 

0.49075 

0.47204 

0.50261 

0.48239 

0.51308 

0.49064 

0.52253 

b-Kja 

0.46607 

0.51308 

0.47972 

0.51792 

0.48607 

0.52253 

0.49115 

0.52694 


Table 5: {M{N) = 124) and . 




0 pi/4 pi/2 3*pi/4 pi 
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Figure 12: 


and Cf" for T e jj ^ 


1 , p = 1 with estimate based on four reference tetrahedrons. 



OL = 

TT 

6 

OL = 

77 

4 

a = 

77 

3 

OL = 

77 

2 

0 

^M,Tr 

i=ir 

C(P 

^M,Tr 

Llp 

C(P 

^M,Tr 

C(P 

^M,Tr 

C(P 

7r/6 

1.09760 

3.78259 

0.96245 

2.71866 

0.91255 

2.27382 

0.93123 

2.05449 

TV/A 

1.09760 

2.43897 

0.96245 

1.78094 

0.91255 

1.50166 

0.93123 

1.38951 

'K jZ 

0.89122 

1.74467 

0.79146 

1.31130 

0.75950 

1.12431 

0.78904 

1.06349 

7r/2 

0.98017 

1.22920 

0.83132 

0.83133 

0.76290 

0.76291 

0.75199 

0.75200 

27r/3 

1.17698 

1.74467 

0.99473 

1.31130 

0.90578 

1.12431 

0.86463 

1.06349 

"i-K jA 

1.35195 

2.43897 

1.14144 

1.78094 

1.03737 

1.50166 

0.98220 

1.38951 


1.65317 

3.78259 

1.39424 

2.71866 

1.26490 

2.27382 

1.19017 

2.05449 


OL = 

TT 

2 

OL = 

277 

3 

OL = 

377 

4 

OL = 

577 

6 

0 

^M,Tr 


^M,Tr 

CjT” 

^M,Tr 

C(P 

^M,Tr 

C(P 

7r/6 

0.93123 

2.05449 

1.07244 

2.39471 

1.21573 

2.95902 

1.47044 

4.21999 

7r/4 

0.93123 

1.38951 

1.07244 

1.64324 

1.21573 

2.01841 

1.47044 

2.80588 

7r/3 

0.78904 

1.06349 

0.91773 

1.27423 

1.04309 

1.50833 

1.26357 

2.11790 

7r/2 

0.75199 

0.75200 

0.86459 

0.86463 

0.98220 

1.12971 

1.19017 

1.67033 

27r/3 

0.86463 

1.06349 

0.96174 

1.27423 

1.08134 

1.50833 

1.30191 

2.11790 

37r/4 

0.98220 

1.38951 

1.07921 

1.64324 

1.20686 

2.01841 

1.44721 

2.80588 


1.19017 

2.05449 

1.29582 

2.39471 

1.44268 

2.95902 

1.72383 

4.21999 


Table 6: {M{N) = 124) and Cj?^ for different 9,aG (0,7r). 


the paper [9] . Constants in functional (embedding) inequalities arise in a posteriori error estimates of the functional 
type (error majorants). The details concerning last application can be found [35l [20l [25l [36l [37l [SU [38] and other 
references cited therein. Below, we deduce an advanced version of an error majorant, which uses constants in 
Poincare-type inequalities for functions with zero mean traces on inter-element boundaries. This is done in order to 
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maximally extend the space of admissible fluxes. However, first, we shall discuss the reasons that invoke Poincare- 
type constants in a posteriori estimates. 

Let u denote the exact solution of an elliptic boundary value problem generated by the pair of conjugate operators 
grad and —div (e.g., the problem dMD-dlll) considered below) and be a function in the energy space satisfying 
the prescribed (Dirichlet) boundary conditions. Typically, the error e := u — v is measured in terms of the energy 
norm ||Ve|| (or some other equivalent norm), whose square is bounded from above by the quantities 


R{v,divq) e dx, / D{Vv,q) ■ V e dx, and / Rr{^{v,q ■ n) e ds, 


Tn 


where H Lipschitz bounded domain, Pat is the Neumann part of the boundary dVl with the outward unit normal 
vector n, and q is an approximation of the dual variable (flux). The terms R, D, and Rtm represent residuals of 
the differential (balance) equation, constitutive (duality) relation, and Neumann boundary condition, respectively. 
Since v and q are known from a numerical solution, fully computable estimates can be obtained if these integrals are 
estimated by the Holder, Friedrichs, and trace inequalities (which involve the corresponding constants). However, 
for H with piecewise smooth (e.g., polynomial) boundaries these constants may be unknown. A way to avoid these 
difficulties is suggested by modifications of the estimates using ideas of domain decomposition. Assume that is a 
polygonal (polyhedral) domain decomposed into a collection of non-overlapping convex polygonal sub-domains 
i.e., 

n:= y nH,»=0, * = 

^ On 

We denote the set of all edges (faces) by Q and the set of all interior faces by Qi^t (i.e., Py G C/int, if by = n Hj). 
Analogously, Gn denotes the set of edges on Pa?. The latter set is decomposed into Pat^ := Pat H dGlk (the number 
of faces that belongs to Gn is ATat). Now, the integrals associated with R and Arw can be replaced by sums of local 
quantities 

/ i?n(u, divg) e dx, and / Rr,^{v,q ■ n) e ds. 

If the residuals satisfy the conditions 


J Rfi^ {v, divg) dx = 0 Vi = 1,..., iV, 


and 


then 


and 


/ Rrfj{v, q-n) ds = 0, Vfc = 1,..., ATat, 


J Rn{v,divq)e dx < C'Qjli?n,(u,divg)||n, ||Ve||n, 


J Rr,^iv,q-n)e ds < C'f;j|i?r„(n,q ■ n)||r^J|Ve||n,. 


(35) 


(36) 


Hence, we can deduce a computable upper bound of the error that contains local constants Cq. and for simple 
subdomains (e.g., triangles or tetrahedrons) instead of the global constants associated with Gl. 

The constant Cq may arise if, e.g., nonconforming approximations are used. For example, if v does not exactly 
satisfy the Dirichlet boundary condition on Pdj, , then in the process of estimation it may be necessary to evaluate 
terms of the type 

J GDiv)eds, /c = 1,..., ATp), 


where Pd^ is a part of Pd associated with a certain Hfc, and Gd(v) is a residual generated by inexact satisfaction of 
the boundary condition. If we impose the requirement that the Dirichlet boundary condition is satisfied in a weak 
sense, i.e., {|Gi:i(n)|}p = 0, then each boundary integral can be estimated as follows: 


I G,,(u)e ds<GP^J|Gz,(n)||p,J|Ve||o,. (37) 

ri:>. 
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After summing up (1351) . (l36l) . and (l37)) . we obtain a product of weighted norms of localized residuals (which are 
known) and ||Ve||o. Since the sum is bounded from below by the squared energy norm, we arrive at computable 
error major ant. 

Now, we discuss elaborately these questions within the paradigm of the following boundary value problem; find 
u such that 


divp + g'^u = /, 

in n. 

(38) 

p = AVm, 

in n. 

(39) 

u = Ud, 

on F£), 

(40) 

AS/u ■ n = F 

on Fat. 

(41) 


Here / G L^(n), F G L?{r n), ud G and A is a symmetric positive definite matrix with bounded coefficients 

satisfying the condition Ai|^p < where Ai is a positive constant independent of The generalized solution 

of (|3^ " (|4T]l exists and is unique in the set Vo + ud, where Vo := |u> G \ u> = 0 on rD|. 

Assume that uGho+Muisa conforming approximation of u. We wish to find a computable majorant of the 
error norm 

lief :=||Ve|li+||gef, (42) 


where || Ve||^ := / AVe ■ Ve da;. First, we note that the integral identity that defines u can be rewritten in the form 


J AVe ■ Vic dx + J f e w dx = J {fw — g^v w — AVu • Vw) da; + J Fw ds, 
Q Q Q r ;V 


Vic G Vn 


(43) 


It is well known (see Section 4.2 in [35]) that this relation yields a computable majorant of |||e|f, if we introduce a 
vector-valued function q G div), such that q ■ n € L^(H), and transform ((43ll by means of integration by parts 
relations. The majorant has the form 


|||e||| < ||£>n(V'c,g)|U-i -f <71^(1;,divg)||n -f C'2||i?r„(i^,g • «)||r,v> (44) 

where Ci and C 2 are positive constants explicitly defined by Ai, the Friedrichs’ constant Cq in inequality 
If ||n < C'nllV^^lln for functions vanishing on F^, and constant in the trace inequality associated with F^r. 
The integrands are defined by the relations 

D{\7v, q) := AVu — q, R{v, divg) := divg + f — g^v, and iipw {u,q • n) := q ■ n — F. 

In general, finding Cq and may not be an easy task. We can exclude C 2 if q additionally satisfies the condition 
q ■ n = F. Then, the last term in (l44l) vanishes. However, this condition is difficult to satisfy, if F is a complicated 
nonlinear function. In order to exclude Ci together with C 2 , we can apply domain decomposition technique and use 
(1351) instead of the global estimate. Then, the estimate will operate with the constants Cq, (whose upper bounds 
are known for convex domains). Moreover, it is shown below that by using the inequalities ([2|) and (l4|), we can 
essentially weaken the assumptions required for the variable q. 

Define the space of vector-valued functions 

H{Q,On,div) := ^q e \ g = G H(Dj, div), {\dWq^ + f - g'^ v\} = 0, \/QiGOn, 

Qj) * ~ ^ ^int; 

{k* • nfe - F|}p„^ = 0, VA; = 1,...,F^}. 

We note that the space F(H, Oq, div) is wider than H(D, div) (so that we have more flexibility in determination of 
optimal reconstruction of numerical fluxes). Indeed, the vector-valued functions in i7(0,div) must have continuous 
normal components on all F^ G Qint and satisfy the Neumann boundary condition in the pointwise sense. The 
functions in F[{Q,Od, div) satisfy much weaker conditions: namely, the normal components are continuous only in 
terms of mean values (integrals) and the Neumann condition must hold in the integral sense only. 

We reform (1331) by means of the integral identity 


/ (g • Vw-I-divg w) da; : 


E 


/ (gi - Qj) 
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which holds for any w G Vq and q G On, div). By setting w = e in (l43l) and applying the Holder inequality, 

we find that 


|||ef < ||i:i(Vn,q)|U-i||Ve|U + divq)||nje - 

+ ^d(9)||e-{|e|}ry Hr, 

r ij G 0int 


Fat^ GFiv 


e — 






where 

rij{q) ■■= life - qj) ■ Iky and pk{q) := \\qk ■ rik - ^llr^v, ■ 

In view of ([1]) and o, we obtain 

lief < ||0(V^,q)fe-x||Ve|U+ ^ ||i?fedivQ)||nfe^J|Ve||n. 

+ ^ ryfeCikj|Ve||n,+ ^ J| Vefey (45) 

Tij-eCint rM'j.GrAr 

The second term in the right hand side is estimated by the quantity ^i(v,q) ||Ve||n, where 

OiSOsj 

We can represent any Gli G On as a sum of simplexes such that each simplex has one edge on dfli. Let fe-feax denote 
the largest constant in the respective Poincare-type inequalities (H]) associated with all edges of SHy Then, the last 
two terms of ((1511 can be estimated by the quantity ^ 2 {v, q) || Ve||n, where 

^l{q)-= (C'^ax)N?. with 77^= jrl{q)+ Y Pkiq)- 

rjjnanj^0 rf^ndQ .^^0 


Then, dlSj) yields the estimate 

III e IIP < |fe(Vu, q)\\A-i IIVe|U + q) + Mq)) II Ve||n, 

which shows that 

|||e||| < ||D(V?;,q)fe-i-b -b 5i2(g)). (46) 

Here, the term 5R2(fj') controls violations of conformity of q (on interior edges) and inexact satisfaction of boundary 
conditions (on edges related to Tat). It is easy to see that K2(g) = 0, if and only if the quantity q-n is continuous on 
Qint and exactly satisfies the boundary condition. Hence, 5R2((?) can be viewed as a measure of the ‘flux nonconfor¬ 
mity’. Other terms have the same meaning as in well-known a posteriori estimates of the functional type, namely, 
the first term measures the violations of the relation q = HVn (cf. (1,491) 1. and 5Ri(u,(7) measures inaccuracy in the 
equilibrium (balance) equation (|38ll . The right-hand side of (1461) contains known functions (approximations v and 
q of the exact solution and exact flux). The constants fe^iax can be easily computed using results of Section [5]111 
Finally, we note that estimates similar to (1461) were derived in [36) for elliptic variational inequalities and in [25) for 
a class of parabolic problems. 
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